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Abstract 

A  gain-scheduling  approach  for  uncertain  Linear  Parameter- Varying  (LPV)  systems 
with  fixed  linear  fractional  relationships  on  a  parameter  set  is  developed.  The  approach 
combines  LPV  theory  based  on  Linear  Matrix  Inequalities  (LMIs)  and  ;u-synthesis  to 
form  a  new  robust  approach  for  large  envelope  control  design.  The  new  approach  is 
used  to  design  an  automatically  gain-scheduled  pitch-rate  controller  for  the  F-16  Vari¬ 
able  Stability  In-Flight  Simulator  Test  Aircraft  (VISTA) .  The  ability  of  the  approach 
to  generate  controllers  for  predicted  Level  1  fiying  qualities  is  illustrated  with  a  high 
fidelity  nonlinear  simulation. 

Keywords:  Aircraft  Control,  Automatic  Gain  Control,  Hqo  control.  Linear  Multivari¬ 
able  Systems,  Structured  Singular  Value 


1  Introduction 


The  challenge  of  most  control  problems  is  to  obtain  as  much  performance  as  possible  in  the 

presence  of  uncertainty  due  to  unmodeled  dynamics  and  plant  variations.  Unfortunately, 

performance  often  must  be  sacrificed  for  robustness  and  vice  versa.  This  typically  happens 
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when  the  physical  plant  varies  significantly  during  operation.  Under  these  conditions,  the 
classical  designer  has  been  forced  to  do  several  point  designs  and  gain-schedule  the  resulting 
set  of  controllers  [1]. 

Several  new  approaches  have  emerged  in  the  last  few  years,  however,  which  may  im¬ 
prove  the  entire  gain-scheduling  process.  These  methods  explicitly  take  into  account  the 
relationship  between  real-time  parameter  variations  and  performance.  Controllers  can  then 
be  designed  for  whole  ranges  of  operating  conditions. 

The  purpose  of  this  paper  is  to  demonstrate  a  new  gain-scheduling  technique  for  un¬ 
certain  LPV  systems.  The  technique  presented  is  a  modified  version  of  the  LPV  method 
developed  by  Apkarian  and  Cabinet  [2,  3].  Their  original  approach  applies  to  systems  of  the 
form 

x{t)  =  A{9{t))x{t)  -f 

y{t)  =  C{0{t))x{t)-\r  D{6{t))u{t) 

where  6{t)  is  a  vector  of  measurable  time- varying  parameters  and  A,B,C  and  D  are  fixed 
linear  fractional  functions  of  6.  In  this  paper,  Apkarian  and  Gahinet’s  LPV  approach  is 
modified  and  mixed  with  techniques  used  in  yu-synthesis  [4,  5]  to  develop  a  robust  synthesis 
method  for  uncertain  LPV  systems  of  the  form 

x{t)  —  A{0[t),5)x{t) B{6[t)^S)u(t) 
y{t)  =  C(9{t),S)x{i)A  D{6(t),5)u{t) 

where  5  represents  the  structured  uncertainty  of  the  system  and  A,  B,  C  and  D  are  fixed 
linear  fractional  functions  of  6  and  5.  The  new  mixed  approach  is  then  used  to  design  a 
pitch-rate  controller  for  the  F-16  VISTA. 

The  remainder  of  this  paper  is  broken  into  3  sections.  Section  2  introduces  some  of  the 
necessary  background  theory  and  the  mixed  LPV//r  approach.  In  Section  3,  the  approach 
is  used  to  design  a  pitch-rate  controller  for  the  F-16  VISTA  and  the  closed-loop  system  is 
tested  with  a  high  fidelity  nonlinear  simulation.  Section  4  summarizes  the  results  presented 
in  this  paper. 
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The  notation  used  throughout  the  paper  is  fairly  standard.  For  real  symmetric  matrices, 

X, X  <0  stands  for  “negative  definite”  and  means  that  all  the  eigenvalues  of  X  are  negative. 
Similarly,  X  >  0  stands  for  “positive  definite”  and  means  that  all  the  eigenvalues  of  X 
are  positive.  The  upper  and  lower  Linear  Fractional  Transformations  (LFTs)  of  a  block 
partitioned  matrix,  M,  with  a  constant  matrix,  N,  are  defined  as 

Fu{M,N)  =  M22  +  M2iN{I-MnN)-'^Mu 
Fi(^M,  =  Mil  T  Mi2N(^I  —  M22X^  ^ M21 

respectively.  For  a  positive  definite  Hermitian  matrix,  L,  is  defined  such  that  F  =  V-F 
satisfies  F*F  =  L.  Finally,  for  an  arbitrary  matrix,  F,  Ker{Y)  represents  the  null  space  of 

Y. 


2  Gain- Scheduling  Uncertain  LPV  Systems 

This  section  briefly  describes  gain-scheduling  for  uncertain  LPV  systems.  The  results  pre¬ 
sented  here  build  upon  the  results  presented  in  [2,  3],  to  which  the  reader  is  referred  for 
further  details  and  proofs.  The  basic  LPV  structure  for  designing  gain-scheduled  Hqo  con¬ 
trollers  is  shown  in  Figure  1. 

The  LPV  plant  is  given  by  the  upper  LFT 


where  P{s)  is  a  known  Linear  Time  Invariant  (LTI)  plant,  0  =  0(t)  is  a  time-varying 
parameter  block  with  the  structure 


0(f)  G  {blockdiag(6>i/,.i, . . . ,  OmFm)}  ■  Oi{t)  G  C  (2.2) 

and  A  represents  a  constant  uncertainty  with  the  structure 

A  G  {blockdiag(Ai, . . . ,  Am)}  :  A,-  G  (2.3) 
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(Note  that  the  Aj’s  associated  with  the  performance  channels  are  purely  fictitious).  It  is 
assumed  the  maximum  singular  value  of  0  and  A  is  less  than  I/7  for  some  positive  scalar, 
7.  The  goal  is  to  find  the  LTI  system,  K,  such  that  the  gain-scheduled  controller  given  by 
the  lower  LFT, 


Fi{K{s),Q)  (2.4) 

ensures  the  closed-loop  system  has  robust  performance. 

Since  the  formulation  presented  below  is  based  on  state-space  matrices,  P{s)  is  repre¬ 
sented  by 

^  Dee  Dei  Db2  \  (  Ce^ 

P{s)=  Die  Dll  Di2  +  Cl  {sI-A)-\Be  Bi  B^)  (2.5) 

D2e  D21  D22  y  y  C'2  y 

where  the  partitioning  is  conformable  to  the  inputs  dg,  d  and  u  and  the  outputs  ee,  e  and  y. 
The  problem  dimensions  are  given  by: 

A  G  IR^^^,Dee  €  M^^^Dn  G  G 

It  is  assumed  that  [A,  B2,  C2)  is  stabilizable  and  detectable,  D22  =  0  and  either  De2  =  0  or 
D2e  =  0.  The  final  assumption  ensures  that  the  gain-scheduled  controller  is  causal  [2]. 

In  [2,  3],  LMI  feasibility  conditions  are  presented  to  construct  an  LFT  controller,  IT,  and 
a  positive  definite  symmetric  scaling,  L,  given  by 

I  Li  L2  \ 

(2.6) 

L3  j 

with  ^20  =  0L2  and  Li^Lz  G  L©,  where 


Lq  =  {Li  >  0  :  LiQ  =  0Li,V0}  C  with  r  =  ^  rj 

i—1 

for  which  the  LTI  system,  Fi(Pa{s),  K{s)),  is  internally  stable  and  satisfies 

(  L^/2  0  \  (  L-1/2  0  \ 

F,(P„(s),/T(s))  <7 

I  0  /  /  0  7/ 
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Equation  2.7  guarantees  Equation  2.4  is  a  7-suboptimal  gain-scheduled  Hqo  controller 
for  the  nominal  LPV  plant.  The  feasibility  conditions  in  [2,  3]  form  a  convex  program  with 
portions  of  the  scaling  matrices  as  the  unknown  matrix  variables.  In  order  to  find  the  scalings 
associated  with  the  minimum  7  in  Equation  2.7,  the  LMI  feasibility  conditions  have  to  be 
solved  iteratively  for  7.  Alternatively,  7  can  be  minimized  subject  to  the  existence  of  the 
real  symmetric  matrices  {R,S)  €  IR"^”,  and  {Lz^Jz)  G  Ta  such  that 


(2.9) 


(2.10) 


5 


^5+7 


(  WIL3WS2  0  \  ^ 

<7 


WlL3Ws2  0 


(2.12) 


0  I3 ) 

\ 

0 

Lz  j 

:  full  column  rank  and 

^  AR  A  RA^ 

RCJ 

Bjs  ^ 

II 

CiR 

-7J3 

biiJz 

Nr 

(2.13) 

y  J3BJ 

mi 

-7-4  ! 

1 

^  A^S  +  SA 

SBi 

C^Lz  ^ 

II 

BjS 

1 

CO 

biiLz 

Nr 

(2.14) 

1 

y  LzCi 

L3DII 

-7L3  j 

N, 


R 


^  Wm  0  ^ 
Wr2  0 
0  I 


Wm 


Ker[  Bj  Dj^ 


(2.15) 


Ns 


^  Wsi  0  ^ 

Ws2  0 
0  I 


Ws2 


=  Ker 


^  C2  D29 


D 


21 


(2.16) 


Bi  =  [  Be  ,  Cl  = 


Ce 

Cl 


i  ^11  — 


Dee  Dei 

Die  Dll 


(2.17) 


^  Ls  0^ 
0  I 


,  J3  = 


h  0 
0  / 


(2.18) 


When  Wr2  and  Ws2  have  full  column  rank,  Equations  2.9  -  2.12  form  a  Generalized 
Eigenvalue  Problem  (GEVP)  (Note:  the  matrix  sums  on  the  left-hand  sides  of  Equations  2.11 
and  2.12  are  independant  of  7  once  simplified)  which  can  be  solved  using  LMI-LAB^^  [6]. 
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If  Wr2  and  Ws2  do  not  have  full  column  rank,  additional  dummy  variables  can  be  used  to 
form  a  slightly  modified  GEVP  [6].  This  modified  GEVP  formulation  was  used  to  solve  the 
example  problem  presented  in  this  paper.  Depending  on  the  dimensions  of  de,  d,  ee  and  e, 
and  the  desired  accuracy  of  7,  however,  it  may  be  more  efficient  to  solve  the  LMI  feasibility 
problem  iteratively  for  ^min-  In  general,  GEVPs  require  more  computational  effort  than  LMI 
feasibility  problems  and  additional  dummy  variables  increase  the  computation  burden. 

Notice  from  Equation  2.7  that  the  channels  associated  with  the  uncertainty.  A,  are  com¬ 
pletely  unsealed.  Adding  scaling  matrix  variables  to  these  channels  ruins  the  convexity  of 
the  above  problem.  The  above  formulation  still  has  value,  however,  since  the  current  ver¬ 
sion  of  /i- Tools  [5]  cannot  handle  repeated  scalar  blocks  or  static  scaling  matrices  efficiently. 
The  deficiencies  in  both  the  LPV  and  //-synthesis  approaches  motivate  a  mixed  LPV  ap¬ 
proach  which  can  exploit  the  structure  of  0  and  A  less  conservatively  while  minimizing  the 
structured  singular  value  of  the  closed-loop  system.  The  following  modified  “D-K”  iteration 
scheme,  dubbed  “D-A'-D”  is  proposed. 


1.  Define  a  scaling  associated  with  the  uncertainty,  La  (5)  G  and  set  it  equal  to 

identity. 


2.  Form  the  inner  D  scaling. 


Di  = 


I  0 
0  La 


(2.19) 


and  the  system 


Ti{s)  =  DiFi{Pa{s),K{s))D;'^ 


(2.20) 


3.  Use  the  GEVP  to  find  an  outer  D  scaling  (the  first  “D”), 

^  Li/2  0 


Do  = 


\ 


0  / 


such  that 


7  = 


DoTiD-^ 


(2.21) 


(2.22) 


7 


is  minimized. 


4.  Design  an  Ho©  controller,  K ,  for  the  system 


p.M 


d2iL~^^^  0 


+ 


\ 


\  ) 


{sI-A)-U  hyL-^!^  b 


using  either  the  DGKF  [7]  central  Hoo  controller  equations  (if  the  plant  is  regular),  or 
the  Hoo  controller  equations  available  in  LMI-LAB^^  [6]. 


5.  Form 


n{s)  =  DoFiiPa{s),K{s))D;'^  (2.23) 

and  over  a  frequency  range  find  scalars,  m§,  and  structured  matrices,  minimizing 


{  mel 

> 

0 

^  mel 

0  ^ 

-1 

IX  = 

To 

\  0 

Ma  ^ 

1 

0 

Ma  j 

6.  Normalize  Ma  at  each  frequency  by  me  and  obtain  La  by  fitting  a  transfer  function 
to  the  normalized  Ma ’s. 

7.  Form  Di  (the  second  “D”)  and  T,-  and  go  to  Step  3.  Stop  when  no  further  reduction 
in  Step  5  can  be  obtained. 

Note  that  the  above  scheme  is  not  guaranteed  to  find  the  global  optimum  solution, 
however,  it  has  worked  well  in  several  practical  examples. 

3  A  Flight  Control  Example 
3.1  The  LPV  Model 

This  section  presents  an  application  of  the  “D-/T-D”  approach  to  the  longitudinal  control 
of  the  F-16  VISTA.  Only  the  longitudinal  motions  of  the  aircraft  are  considered.  Roll  and 
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yaw  angular  velocities  and  the  sideslip  angle  of  the  aircraft  are  assumed  to  be  zero.  The 
standard  short-period  equations  of  motion  [1], 


M, 


«  / 


(a\ 

9  / 


(3.1) 


are  used  for  design  purposes.  In  Equation  3.1  a  represents  the  angle-of  attack,  q  is  the 
pitch-rate,  and  5^  is  the  elevator  deflection.  The  effects  of  the  trailing  edge  flaps  on  the 
longitudinal  dynamics  are  ignored.  At  trimmed  level  flight  the  dimensional  coefficients,  Zai 
Ma,  Mg,  Zs^,  and  Ms^  depend  mainly  on  Mach  number,  M,  and  altitude,  h. 

The  Flight  Dynamics  Directorate  of  Wright  Laboratory  has  a  high  fidelity,  six  degree 
of  freedom,  nonlinear  simulation  model  for  the  F-16  VISTA  [8].  This  simulation  contains 
accurate  models  of  the  propulsion  system,  actuators,  sensors,  disturbances,  loadings,  atmo¬ 
sphere,  equations  of  motion,  sensors,  and  an  aerodynamic  database  for  a  large  range  of  Mach 
numbers,  altitudes,  and  angles  of  attack.  This  nonlinear  model  was  trimmed  and  linearized 
at 


M 


0.35  0.45  0.55  0.65  0.75  0.85 


h 


1000  ft  5000  ft  15000  ft  25000  ft 


(3.2) 


The  data  was  then  fit  with  polynomial  functions  of  M  and  h.  The  polynomial  expressions 
for  the  coefficients  are: 


=  (2.20e  -  1)  -  (4.10e  -  7)h  -  (2.60e  +  0)M -f- (5.15e  -  5)M/i  (3.3) 

Ma  =  (1.71e  +  l)-(8.07e-4)h-(6.84e-l-l)M-t-(3.31e-3)Mh 

-h(5.62e  -h  1)M^  -  (2.92e  -  3)M^h  (3.4) 

Mg  =  -(2.28e  -  1)  4- (7.06e  -  6)h  -  (2.12e -h  0)M -h  (4.86e  -  5)Mh  (3.5) 

Zs,  =  -(1.38e-3)-|-(8.75e-8)/i-(3.40e-l)M-t-(7.98e-6)Mh  (3.6) 

Ms,  -  -(8.16e  +  0)  +  (1.73e-4)/i-l-(4.06e-l-l)M-(8.96e-4)Mh 

-(9.93e  +  1)M^  +  (2.42e  -  3)M^h  (3.7) 
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The  design  region  includes  the  flight  conditions, 


Me 


0.4  0.8 


h  e 


mo  ft  25000  ft 


(3.8) 


which  covers  roughly  one  third  of  the  VISTA’s  Ig  sustained  maneuver  envelope  at  military 
power.  Since  the  LPV  approach  presented  above  is  based  upon  the  Small  Gain  Theorem  [9] 
it’s  convenient  to  express  Equations  3.3-3. 7  in  terms  of  the  normalized  variables 


M-0.6 


-1  1 


„  h-  15000  ft 
Sh  =  . : . 6 


-1  1 


10000  ft 

Letting  Om  and  Oh  equal  the  order  of  the  mach  and  altitude  fits,  respectively,  each  normal¬ 
ized  polynomial  expression  can  be  converted  to  an  LET  expression  in  terms  of 

0  =  blockdiag{SMInM-,  ^hlnh)  (3.9) 

where  nM  —  OM{Oh  +  1)  and  nh  =  Oh  by  writing  the  dimensional  coefficients  as 

Ok  =  Q22  +  ^21(0-'  -  Qn)-^  =  E„(g,  0)  (3.10) 

The  partitioned  matrix  Q  is  equivalent  to  a  realization  of  a  multi-dimensional  system,  so 
its  “order”  (the  size  of  0)  can  be  reduced  without  penalty  by  removing  “uncontrollable 
and  unobservable  states”  associated  with  5M  and  5h  in  turn.  While  this  approach  doesn’t 
guarantee  the  resulting  0  will  have  the  minimum  possible  size,  it’s  easy  to  implement,  and 
works  well  in  practice. 

3.2  Problem  Setup 

The  problem  is  to  design  a  pitch-rate  controller  for  the  F-16  VISTA  which  provides  robust 
command  tracking  with  predicted  Level  1  handling  qualities  [10].  Time  domain  specifications 
for  the  pitch-rate  response  are  shown  in  Figure  2  and  Table  1.  Note  that  the  “rise-time” 
parameter,  St,  varies  with  the  true  velocity  of  the  aircraft  (in  ft/sec).  In  order  to  achieve  Level 
1  flying  qualities  the  model  matching  control  structure  in  Figure  3  is  used.  The  reference 
model  is  the  second  order  system, 

42 


s2  +  2(0.5)(4)s+42 


(3.11) 
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which  meets  Level  1  flying  qualities  over  the  entire  design  envelope.  In  keeping  with  the 
baseline  VISTA  control  design,  one  reference  model  was  chosen  so  that  speed  of  the  pitch-rate 
response  is  virtually  the  same  over  the  entire  design  envelope.  The  maximum  allowable  pitch- 
rate  command  (before  the  actuators  saturate)  varies  with  flight  condition.  Alternatively,  a 
parameter-varying  reference  model  could  be  used  in  the  design. 

The  design  plant  in  Figure  3  also  contains  a  first  order  approximation  of  the  actuator 
dynamics. 


20.2 
5  -f-  20.2 


(3.12) 


input  multiplicative  uncertainty,  sensor  disturbances,  and  weights  on  the  elevator  deflection 
rate,  uncertainty,  noise,  and  pitch-rate  error.  The  various  design  weightings  are  shown  in 
Figure  4.  The  control  weighting,  IFc,  reflects  the  maximum  allowable  elevator  deflection  rate 
of  70  deg/sec  and  the  noise  weighting,  Wn,  depicts  a  5%  error  expected  in  the  measurements. 
The  performance  and  uncertainty  weightings,  Wp  and  are  chosen  as  lowpass  and  highpass 
filters  respectively  to  provide  the  standard  performance/robustness  tradeoff.  The  exact 
values  for  Wp  and  Wu  were  selected  after  several  design  iterations. 

A  plot  of  the  maximum  singular  value  of  the  closed-loop  system  after  Step  5  in  the  7 
step  “D-K-D”  process  for  the  and  4*^  iterations  is  shown  in  Figure  5.  The  maximum 
scaled  singular  value  after  the  1®*  iteration  (corresponding  to  the  unmodified  LPV  approach) 
is  1.48.  After  4  iterations  the  upper  bound  to  [j,  converges  to  0.973.  The  13*^^  order  controller 
found  after  the  4*^  iteration  was  reduced  to  7’*^  order  by  truncating  the  fastest  modes  and 
residualizing  the  resulting  system. 


3.3  Nonlinear  Simulation 

A  high-fidelity  nonlinear  parameter-varying  step  response  of  the  system  is  shown  in  Figure  6. 
The  aircraft  is  initially  perturbed  from  trimmed,  level  flight  at  M  =  0.6,  and  h  =  15000  ft. 
The  pitch-rate  response  demonstrates  predicted  Level  1  flying  qualities. 
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4  Summary 


A  new  LPV/;U  design  approach  was  presented  and  used  to  design  an  automatically  gain- 
scheduled  pitch-rate  controller  for  the  F-16  VISTA.  The  design  was  tested  with  a  high  fidelity 
nonlinear  simulation  and  Level  1  flying  qualities  were  predicted  from  the  time  responses.  The 
method  presented  takes  advantage  of  familiar  concepts  in  linear  control  theory  (i.e.  LFTs, 
and  structured  singular  value  synthesis)  and  new  convex  optimization  tools  (LMIs).  The 
major  disadvantage  in  the  approach  is  the  lack  of  direct  way  to  incorporate  realistic  rate 
limitations  on  the  changing  parameters  of  the  system.  Overall,  however,  the  new  approach 
is  promising  and  warrants  further  research. 
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Table  1:  Pitch-rate  handling  qualities  specifications 


Parameter 

Level  I 

Level  II 

ti  max 

0.12  sec 

0.17  sec 

Agz/Aqi 

0.30 

0.60 

At  max 

SOO/Vt  sec 

1600 /Vt  sec 

At  min 

9/Vt  sec 

3.2/Vr  sec 
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